Hamiltonian and Linear-Space Structure for Damped Oscillators: I. General Theory 
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The phase space of N damped linear oscillators is endowed with a bilinear map under which 
the evolution operator is symmetric. This analog of self-adjointness allows properties familiar from 
conservative systems to be recovered, e.g., eigenvectors are "orthogonal" under the bilinear map 
and obey sum rules, initial- value problems are readily solved and perturbation theory applies to the 
complex eigenvalues. These concepts are conveniently represented in a biorthogonal basis. 
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I. INTRODUCTION 

Conservative systems are described by self-adjoint 
time-evolution operators "H on a linear space, with well- 
developed mathematical tools. Eigenvectors are orthogo- 
nal and complete, and eigenvector expansions are unique, 
with coefficients given by projection. Expansion of the 
initial data solves the dynamics, with each term simply 
acquiring a phase e _lWjt , where u)j are the eigenvalues. 

Indeed much of our very language, e.g., normal 
modes, energy levels, transition frequencies, photons and 
phonons, is rooted in eigenvector expansions. These ex- 
amples already suggest that the present ostensibly clas- 
sical problems are also relevant in the quantum domain. 

However, these concepts are ruined by dissipation, 
since TL is no longer self-adjoint. This paper shows, for 
ohmically damped oscillators, that many properties well 
known from the conservative case can be resurrected, in 
terms of not the standard inner product, but a bilinear 
map under which TL turns out to be symmetric. This 
symmetry is analogous to self-adjointness, and leads to 
the orthogonality of eigenvectors under the bilinear map, 
the usual ready solution of initial-value problems and 
Rayleigh-Schrodinger perturbation theory (RSPT) in the 
familiar form but applicable to complex eigenvalues. 

Serious microscopic discussion of dissipation began 
with a single linear oscillator coupled to an infinite bath, 
which is then eliminated from the equations of motion or 
path integral [lj . A huge body of works generalize this 
to many nonlinear coupled oscillators @- Here, however, 
we first deal with N coupled linear, classical oscillators, 
with coordinates 4>(a), a — 1, . . . , N, described by 

(%4>{a) + Y(a, L3)dt<t>{P) + K(a, p)<j>(fi) = (1.1) 

(summation convention for Greek indices), with Y a sym- 
metric non-negative damping matrix and K a symmet- 
ric positive matrix of force constants One can ob- 
tain (|l.lfl by elimination of a bath, giving the dynamics 
the necessary quantum pedigree; but until Section [VII it 
can simply be regarded as phenomenological. The non- 
negativity of r may be relaxed, e.g., in continuum models 
of the electromagnetic field in gain media. 

Section[H]sets up the linear-space formalism. However, 
eigenvector expansions would not be useful unless there 



is a convenient projection; for this purpose, the key bilin- 
ear map (tp, <p) between state vectors in 27V-dimensional 
phase space (leading to a biorthogonal basis) is intro- 
duced; examples are given in Section lHfl The application 
to time-independent perturbation theory (Sect ion llV|) for 
the complex eigenvalues and eigenvectors proceeds by 
mimicking RSPT for conservative systems. Section 
briefly discusses the continuum limit N — > oo, recover- 
ing results for waves in open systems, which have been 
studied extensively 0. 

Section IVII sketches the application to thermal and 
quantum physics, and to nonlinear oscillators. As a fur- 
ther outlook, Appendix [S] considers the separation of 
time scales if one oscillator is light. Adiabatic elimination 
of its momentum yields effective constrained dynamics, 
giving an example where the methods still apply also 
when phase space is odd-dimensional. An extensive ac- 
count can be found in Ref. p|. 

The formalism assumes that (a) although TL is not self- 
adjoint, its eigenvectors fj are complete, and (b) for all j, 
(fj,fj) 7^ 0. These are in fact related, and violated only 
at critical points (with measure zero in parameter space) 
where eigenvectors merge (Section HTHl and Ref. @). A 
brief treatment focusing on "excess noise" in open sys- 
tems can be found in Ref. 0. 



II. FORMALISM 

A. State space and evolution operator 

We cast (|l.lf) in first-order form by introducing the 
momenta <fi(a). Denote the state in phase space as 



2iV 



= 



,2N 



) T = 5> Mj , (2.1) 



where Uj are the basis vectors of this canonical represen- 
tation, and boldface denotes 2A r -vectors. Then 



dt4> = —iTL<p , 
/ 



H . = i 



-K -r 



(2.2) 
(2.3) 
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with I the N x N identity. The dots specify whether the 
indices are contravariant or covariant (cf. Section 111 D|) . 
Eigenvalues are the roots of 

J(w) = det(W. - u) , (2.4) 

and right eigenvectors fj are defined by 

Hfj = LJjfj . (2.5) 

Left eigenvectors /■? will be introduced later. For fj, 
(|2.2|) implies the time-dependence e~ tuJjt , so 



-IWjfj 



(2.6) 



and one can simply solve l|l.l|) in the frequency domain 
for the coordinates, without reference to the momenta: 

[-^S(a, (3) - iLo 3 T(a, (5) + K(a, 0)] fj (0) = . (2.7) 



With dissipation, ImUj < [8|. Conjugating 1)2. 5)1 . 
it follows that u)_j = —uj* is also an eigenvalue, with 
f_j oc /* (see below for normalization and phase). Thus, 
eigenvalues are paired, except for so-called zero-modes 
with Re ujj = [3 ■ Indeed, the determinant of [• • ■ ] in 
(|2.7() has 2 N roots Wj ; if T ^ then typically also all w| 
are different, showing that an eigenvector expansion can 
only be sought in phase, not coordinate, space. 



B. Inner product and bilinear map 



In phase space, the standard inner product is ^(j 



(-010) = 0(a)*0(a) + ip(a)*ct>(a) , 



(2.8) 



with (0|0> > 0. Unfortunately, (V|W0) ^ (cf>\Hip)* , so 
(fk\fj) 7^ even if uj k ^ io r 

Our key concept is the bilinear map: 

(0, 0) = i fy(a)4>(a) + $(a)0(a) + 0(a)r(a, /3)0(/3)] 
= (0,0). (2.9) 



The diagonal entries (0,0) are not positive definite, not 
even necessarily real [Tl|. By (|2.6f) . (fj,fk) can be ex- 
pressed in terms of the coordinates alone. 

We will justify (|2.9|) by the ensuing properties, cru- 
cially 1)2.10)1 below; still, let us motivate it immediately. 
First consider the conservative case T = 0. A structure 
ijxj) + " t P4 > is n °t dimensionally correct |10| , but 0</> + 00 
is. Moreover, with 1)2. 9fl an eigenvector pair with coordi- 
nates fj{a,t) = / J (a)e~ 4 ^* and f-j(a,t) = / i (a)*e 1 ^* 
will have opposite momenta, hence a zero bilinear map. 
This would not hold if one vector is conjugated. When 
there is damping, the last term 0r0 is the only other 
type with the correct dimension; with a unit coefficient 
as in 1)2. 9fl . it ensures 



(0,H0) = (0,H0) 



(2.10) 



In the proof, r in TC cancels T in (|2.9|l . This analog of 
sclf-adjointness, as usual, yields "orthogonality" 



(/i,/fc) = if cjj^u k , 



(2.11) 



including pairs uj-j, Wj- Taking suitable linear combi- 
nations, (|2.11|) can be extended to level crossings, i.e., 
multiple independent eigenvectors at one eigenvalue. 



C. Normalization and phase 

If (fj,fj) 7^ (cf. Section ITl Hp . completeness implies 

(2.12) 



This normalization-independent form is appropriate if 
small denominators must be handled for (fj,fj) ~ 0. 
Also, the usual convention for the conservative limit is 
fj(a)fj(a) — 1, i.e., (fj,fj) = 2uij [incidentally explain- 
ing the conventional factor i in 1)2.9)1 ]. Thus, connection 
with this limit is best expressed via 1)2.12)1 . 

However, in most other circumstances, it is more con- 
venient to adopt the normalization (and phase) conven- 
tion (fj, fj) = l for all j, in which case 



(/i,/*) = * 



jk 



(2.13) 



Unless otherwise specified, expressions that are not man- 
ifestly normalization- independent conform to 1)2.13)1 . 
Conjugate eigenvectors have the phase relationship 



f—j = ±ifj 
in order to satisfy flTC^ pi ]. 



(2.14) 



D. Metric 



In the u-basis [l3|, (0,0) = Yli j 0*5ij '0 J m terms of 



9ij = (ui,Uj) 



This, and the inverse metric 



.9 = -« 



r i 

I 



I 

1 -I 



(2.15) 



(2.16) 



are used to lower/raise indices as usual. This ex- 
plicit form shows that g.. is non-singular. Moreover, if 
(0, 0) = for every 0, then = (g" 0* , 0) = J2j \^\ 2 ' 
implying = 0. In the eigenbasis, gij = (/,, fj) = 5ij. 

The crucial symmetry 1)2.10)1 translates into the sym- 
metry of TL.. = g.. TL' . . Explicitly, 



H.. = 



K 
—1 



(2.17) 
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in the evaluation of which [cf. 1(2.3(1 and 1(2.15(1 ] the can- 
cellation of r becomes apparent. Incidentally, Ti... can 
be expressed directly in terms of the matrix elements 
{tp,Hcf)) = ^Kijfl , and using 12.17( 1 gives 

(0, H4>) = -0(a)0(a) + <f>(a)K (a, 0)<f>(p) 

= -2L[<f>] , (2.18) 

where L can be interpreted (at least in the non-dissipative 
limit) as the Lagrangian. 

Although r disappears from (|2.17|) , physical quantities 
of course do depend on the damping. For example, u>j 
is the stationary value of (</>, 7i0)/(</>, </>) attained when 
4> = fj. The numerator 1(2.18(1 is T-independent; but the 
denominator carries the nontrivial dependence on T. 

E. Duality 

For any basis {vj}, there is a dual {«•'} such that 
(vi\vk) — 5 J \ 14]- For the canonical basis, u n = u n , 
so that the 2N x 2N coefficient matrices (in the u-basis) 
(vj) n and (v-*)^ are each other's inverse. However, for 
the eigenbasis one can bypass this inversion, since ((2.13|) 
and 1(2.15(1 give (p) m = E n [9mn(/i)™]*, that is, 

f =Vf j = \g..f j ]* . (2.19) 

This is extended to any vector <p, i.e., 2?0 = [.<?•• 0]*- 

The duals are left eigenvectors 11) fi — u*fi. In terms 
of these, = Yl j fj (f J 10) and, in an obvious shorthand, 

/ £/j /' • • (2-20) 

3 

with I the 2N x 2N identity. The dual basis leads to a 
consistent use of contravariant and covariant indices. 

F. Sum rules 
Writing out 1(2.20(1 in components yields four sum rules: 



are paired, obeying 1(2.14(1 . Thus, e.g., f-j(a)f-j(j3) — > 
—fj(a)fj(/3), making ((2.21a(l trivial. The remaining 
((2.21b() reduces to the familiar sum rule for conservative 
systems (cf. Section lH Cl for the different normalizations). 



G. Time evolution 
In terms of l(2.12|l . the dynamics is formally solved as 

0(i) = V^fj > ^ = (fiMt=0)) ■ (2-22) 

i 

This familiar construction would not be possible with the 
product ((2.8(1 . which would have required inverting the 
2N x 2N matrix (fj\fk) — nontrivial if N is large. Of 
course, to use ((2.22(1 . one must know the fj, but when 
only a few / 3 's dominate the dynamics, the burden of 
finding those (e.g., iteratively) is comparatively small. 

In operator form [cf. 1(2.20(1 ], the Green's function is 
4>{t)6(t) — G(t)4>(t=0) (6 is the unit step function), or 

(H - w) G{uj) = -iX . (2.23) 

With the eigenvectors complete, explicitly one has 

«=(|p2 I-)=p7^ (/;>•), (2-24) 

Fourier inversion of which reproduces 1(2.22(1 . 

Using ((2.3(1 . in the canonical basis we can also write 

involving the inversion of only an N x N matrix, in con- 
trast to 1(2.231) . The first factor in ((2.25(1 corresponds to 
the differential operator in l(2.7|) . 






3 


(2.21a) 


I 


= ^Z^jfj®fj > 
3 


(2.21b) 





3 


(2.21c) 





= I> J Vi®/i + ^i/i®(r/ J 0] • 


(2.21d) 



Here, a®b stands for the matrix with elements a(a)b((3). 
We have verified all these in examples. 

For r -> 0, l|2.2lajl . i|2.2lc|) and l(2.21d(l become vac- 
uous. Namely, all Rewj 7^ for small T, so the fj 



H. Exceptions 

Our formalism relies on (a) the completeness of eigen- 
vectors, and (b) (fj, fj) 7^ for all j. These assumptions 
are related. Imagine that system parameters are tuned, 
and the 2N eigenvectors start off (e.g., for T = 0) com- 
plete. Then (a) is violated when two (or more) eigenvec- 
tors merge, e.g. f k -> fj, so (fj,fj) = lim (/_,-,/*,) = 0. 
Conversely, if (b) is violated, fj is orthogonal to every 
eigenvector; if these would be complete, then fj is orthog- 
onal to all 0, contradicting fj ^ 0, cf. below ((2.1tj|) 15]. 

These conditions are violated only for a zero-measure 
set of parameters. Section IlIII contains examples; the ex- 
tension to these critical points will be given separately [|| . 
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III. EXAMPLES 



A. Single oscillator 



(u)j + 2i"fUj)fj(a). But K has a complete real eigensys- 
tem {Qpfj}, with fij the frequencies for 7 = 0. Now 



Although the theory is most valuable when N 1 
(otherwise brute force suffices), it is already nontrivial 
for N = 1. Let K = k and T = 27. The eigenvalue 
equation J(w) = [cf. 12.4|l ] leads to u>± = ±J7 — 17, 
where SI = \/k — 7 2 , and 



1 



The bilinear maps are (/+, /-) = and 
(/±,/±)=±2c|fi. 



(3.1) 



(3.2) 



At the critical point k = fc* = 7 2 , (a) the eigenvalues 
merge: w + = cj_ = —17; (b) from l|3.1|l the eigenvectors 
merge, leaving only one in the 2-d space; and (c) from 
(|3.2() the diagonal bilinear map vanishes. This example 
also explains why eigenvector merging is called criticality. 



Two oscillators 



Let N ~ 2 with Q 



K 



4 
-2 



4 7 
2 7 



(3.3) 



As 7 increases, one mode pair becomes overdamped for 
7 > 7*1 = 0.8599, and the other for 7 > 7*2 = 2.1031. 
Except at these critical points, the eigenvectors are com- 
plete, and all bilinear maps etc. can be found explicitly. 
A perturbation of (|3.3|) will be studied in Section liV Al 



C. Chain of masses 

Next consider a chain of unit masses with displace- 
ments 0(a) from the rest positions x a = aA. Mass a is 
tied to x a and its nearest neighbors by springs of force 
constant V(a) and k respectively, and is moreover sub- 
ject to a damping force — r(a)<9 t 0(a). That is, 

(2k + V(a) if \a-/3\ =0 
K(a,l3) = I -k \{\a-(3\ = l , 

[0 if jor — /3| > 1 



T( a ,/3)=T(a)S(a,/3) 



(3.4) 



In the obvious continuum limit A — > 0, fcA 2 = 1, the 
system obeys the (generalized) Klein-Gordon equation 



[d 2 + T(x)d t -d 2 x + V{x)]cj>{x, t) = , 



(3.5) 



discussed in Section IV Al 

Let < a < N+l, 0(0) = </>(N+l) = 0. In the spe- 
cial solvable case T(a) = 27, l|2.7|l reads K(a, f3)fj(fl) = 



-«7 



(3.6) 



while the eigenvectors remain the same. Thus this model 
is in effect a superposition of many oscillators, each as 
in Section fill Al For any 7, at most one of these can be 
critical. This example is simple because [K, T] = 0. 

Let us comment briefly on the doubling of modes ^(| . 
For r = 0, the Vt 2 are real; but each of these splits into 
two complex w 2 when r ^ 0. Modes are doubled because 
expansions now refer to 4> and simultaneously. 

The /j-basis for the [imjdamped system is orthogonal 
under 12.911 [ l|2.8|l ]. This is consistent, since 



±2 V / 3 2 -f/ 3 W/ J («) 



(3.7) 



equals the conventional norm for the conservative system 
except for the prefactor, vanishing at critical damping. 

To be even more specific, let V(a) — 0, so that the un- 
damped case is essentially the familiar model for acous- 
tical phonons, with fi 2 = 2fc{l — cos[j7r/(A-t-l)]} and 
fj(a) = Cj sin\jira/(N+l)] for 1 < j < N. Setting |3~7jl 
to unity, one finds cj = [(ujj+i^lN+l)]^ 1 / 2 . 



D. Inhomogeneous damping 

A non-trivial solvable model with [K, T] 7^ can be 
constructed as follows. Let the masses be a = 0, . . . , N, 
with again 0(0) = but with <j)(N) free, so there are 
N dynamical variables, with force constants as in H3.4(l 
except that 



K(N,N) = k + V(N) . 
Let only the final mass be damped: 

T(a,(3)= 1 S(a,N)5((3,N) 



(3.8) 



(3.9) 



We only consider V = 0, returning to finite V in Sec- 
tion lV Bl For N — > 00, it is then straightforward to find 



Nuij J — i artanh(7/v / fc) , 7 < \fk : 

yk M7T — i arcoth(7/vfc) , 7 > \fk . 



(3.10) 



The divergence as 7 — > \fk occurs only for iV — * 00. For 
finite N ^> 1, (|3.10(l fails very near 7 = Vk, and the true 
u>j are continuous in 7. In particular, for 7 = yk one has 
uj = {Vk/N)[(j+ X )n - iilniV], where X = |, |, and \ 
for -N < j < - In A, \j\ < In A, and In A < j < A 
respectively. 
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IV. PERTURBATION THEORY 

The familiar RSPT can be transcribed, everywhere re- 
placing the usual inner product with the bilinear map 
(|2.9|) . However, this assumes that the bilinear map itself 
is unperturbed, so we only consider, for small real e, 



H = H + eAH 



AH = i 




-AK 



(4.1) 



A. Simple modes 

The changes in the eigenvectors up to first order, and 
in the eigenvalues up to second order, are 01 



Af j = eJ2n 



10 j — LUk 



Au; j = e(f j ,AHf j 



E 



(fk,AHf 3 ) 

UJj - UJk 



(4.2) 
(4.3) 



In the e 2 -term, the matrix elements are squared without 
taking absolute values. Higher orders are easily writ- 
ten down. Here and below, {u>j,fj} is the unperturbed 
eigensystem. In view of (|4.1() . one simply has 

{fk, KHfi) = f k (a)AK(a, 0)fj(0) = (AK) kj . (4.4) 

Consider the example in Section 1111 Bl with 7 = 0.5 
for H , and set K(l, 1) h-> K(l, 1) + e. Figure □ shows 
Aui(e) for the least damped mode u>\ = 2.267 — 0.699i. 
The points are roots of J(w); the line is given by l|4.3|l . 
the intercept (slope) verifying its first (second) term. The 
remaining error (a quadratic in the plot) reveals an e 3 
term in Alo. RSPT now gives Im Alo correctly as well. 

Re-write the first-order shift allowing arbitrary nor- 
malization of fj, as A(ujj)/c = 2ujj(AK)jj/(fj, fj). 
Agreement with the conservative limit is now apparent. 
This also shows that the bilinear map is no mathemat- 
ical artifact, since it relates directly to the frequency 
shift — (AK)jj being the expected expression. Inciden- 
tally, there could be large shifts if fj) is small; the 
interesting consequences in optics [lg are discussed in 
detail elsewhere 0, • 

Since ImAujj — Xm<^ ^ e ") * s nove l m the 

generalized RSPT, one further aspect deserves mention. 
Suppose the state j is weakly damped with quality fac- 
tor <2, but other states are strongly damped: "fj/10 ~ 
Q -1 <C 1 and Y/oj ~ 1; here jj = |Imwj|, 7' is the typ- 
ical llmwfel for other fc, and w is a typical |Reu;;|. One 
then has = lm[e(AK)jj] ~ ejj [because fj € R if 
7, -» 0], but 5P) ~ Im[e 2 (AX)2 3 /(^-^)] ~ e 2 7 '. The 

ratio 5^/5^ between energy loss through transitions to 
other states and direct change to the damping thus is 
~ e<2, which could be larger than one if « £ « 1 
(very narrow resonance). 
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FIG. 1: The real (a) and imaginary (b) parts of the eigenvalue 
shift Alo for the system in 1)3. 3|l (with 7 = 0.5), due to a 
perturbation adding e to K(l, 1). Straight lines are second- 
order RSPT; points are roots of the characteristic polynomial. 



B. Perturbation around level crossings 

Suppose, as parameters are tuned, M eigenvalues 
merge. Generically, the corresponding eigenvectors also 
merge (criticality) , so the remaining ones are incom- 
plete 0, El- Instead, here we consider level crossing, 
i.e., at an Afth-order zero of J(to) there are M indepen- 
dent eigenvectors. Again exactly as in the conservative 
case, the (typically small) M xM block is explicitly diag- 
onalized, leaving a perturbation without intra-block el- 
ements. Thus zero denominators lo^ — t-Oj are avoided, 
and the formulas in Section llV Al apply: we omit the de- 
tails. The same method avoids small denominators near 
a level crossing, where a group of eigenvalues (but not 
the eigenvectors) are close. 

Near-crossing leads to level repulsion and attrac- 
tion [2(J, which can now be stated precisely. Focus on 
near-degenerate states j = 1,2 and remove the aver- 
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age eigenvalue. In the associated subspace, (Ho) 
diag(<5, — S). Perturb this as in (|4.1|) . with 



AH' . 



k 
k 



(4.5) 



where k — (AK) 12 as in (|4.4|l . Then the eigenvalues are 



CO 



±^S 2 + (ek) 2 



(4.6) 



2+A8 
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K= 


2 




r = o 


i 


• 


V o 


2 












(\6\, \e\ <C 1, but e/S is arbitrary). Note that the phys- 
ically relevant quantity is k 2 , meshing with the normal- 
ization (|2.13|) . which fixes A: up to a sign. 

// the two modes correspond to two near-identical os- 
cillators (as is always the case for N — 2), then / 2 2 and 
hence k are real; cf. Section IIII Al One gets the oft- 
quoted [2(j result that frequencies repel but widths at- 
tract. This also holds for two modes which are both near 
crossing and weakly damped, by l|2.13[) and (|4.4|l 21] . 

However, the general case is not as simple, because K 
and r need not commute. For instance, take 



(4.7) 



6 lifts the degeneracy at to = Perturbation by 

^12 Ml3\ 

AK = | /ii2 (4.8) 
V M13 J 

[with real and O(l)] yields ATL as in 14. 511 . where 

can have any phase. Suf- 
ficiently near the level crossing — the region of interest — 
k can be taken constant as e runs through the reals. Sim- 
ple geometry then shows that the locus of to in 14.6fl is 
part (by the restriction e 2 > 0) of an orthogonal hyper- 
bola, with (un)stable direction ik (k). 



V. CONTINUUM MODELS 

The continuum limit N — > oo is relevant to classical 
and also quantum fSection lVl|l fields, e.g., electrodynam- 
ics in an absorbing medium. Although the examples be- 
low are 1-d, the theory also applies to three dimensions. 



A. String subject to viscosity 



Consider the example in Section Till CI In the obvious 
continuum limit A — ► 0, N — * oo with a = (JV+1)A 
fixed, and without loss of generality scaling kA 2 = 1, 
the dynamics obeys (|3.5() for < x < a, with </>(0) = 
4>(a) — 0. This models a string subject to viscosity. (Gen- 
eralization to0<x<ooor— oo<:r<oo, in both cases 
with 4> vanishing at both ends, is straightforward.) 



In terms of the momentum 4>(x) = dt<j){x) and the two- 
component field 4> = ((/>, <^) T , one has id t <fi = Ti.(p, with 

H W= l (d 2 -V(x) -T(x)J ■ (5 - 1} 
The bilinear map is (up to an irrelevant overall factor A) 



(ip, <p) = i (ip<j> + i)(j> + ipT<f>) dx 



Again, H is symmetric under (|5.2|l . 

The eigenvalue equation is Ttfj = Ujfj, or 

[-to] - tc^rOr) - dl + V(x)] fj(x) = . 



(5.2) 



(5.3) 



This dissipative Sturm-Liouville (DSL) problem could be 
generalized further by — d 2 i— > —d x p(x)d x for suitable p. 
In the familiar case T = 0, each V gives a complete set 
{fj} v , orthogonal under the conventional inner product. 
Now, each (V,F) gives a set {fj} vr , orthogonal under 
(|5.2() and expected to be complete, with the convergence 
of infinite two-component sums guaranteed if (V, T) are 
suitably smooth and bounded. Such DSL systems are in- 
teresting both for, e.g., functional-analytic aspects such 
as completeness and convergence, and for the specific or- 
thogonal functions that emerge. 

If T(x) = 2 7 , (E3J) becomes (V - d 2 x )jj = fl 2 f j: where 
Qj and LOj are related by (|3.t)|) . and V — d 2 has a complete 
real eigensystem {il 2 , fj}. The norm (fj, fj) equals the 

conventional one Jq/J dx up to a pre-factor ±2V^ 2 — 7 2 
[cf. (|3.7|) ]. which vanishes at critical damping. 
All these parallel the discrete version in Section UlI CI 



B. Waves in open systems 

Next consider a similar string on < x < oo, seg- 
mented into the "system" x < a and a trivial \V = T = 0] 
"bath" x > a, into which the waves are outgoing: 



{d t + d x )(j){x>a) = 



(5.4) 



which follows automatically if the bath is at rest initially. 
So the "system" is to be solved with the same outgoing 
(rather than nodal) condition at x = a. For simplic- 
ity, also take r(cc<a) = so that damping occurs only 
through the escape of waves, not through the field equa- 
tion. This models a laser cavity with output coupling. 

To see how also this model is a limit of 111.11 . set 
x = Aa, with A = a/N, < a < N, and ^(0) = 0. 
For a < N, the discrete dynamics can be written down 
immediately. For the last particle, reference to the next 
one just outside the system can be eliminated by l|5.4l) as 
[(f)(N+l) - cj)(N)]/A = -d t cj)(N), resulting in 



d 2 t + x + V(N) 



<t>(N) 



^(jV-l)-^(iV) 
A 2 



(5.5) 
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These exactly correspond to (|3.8|) - (|3.9[1 for k — A~ 2 and 
7 = A -1 — without incoming waves, the outside merely 
acts as a source of ohmic damping on the right endpoint. 
The continuum limit of the bilinear map reads 



00) 0) = i 



(00 + ip<t>) dx + ip(a)(f>(a) 



(5.6) 



(again up to a factor A), i.e., (|5.2[l with T(x) — 5(x—a). 
The vital symmetry of H. under l|5.6l) follows as a limit 
of H2.10|) . For a direct proof integrate the — C* 2 , term by 
parts' the resulting surface term cancels the one in (|5.6|) 
All results then carry through as before. 
Our extensive study of such open wave systems mo- 
tivates much of the present work (notably the ansatz for 
the bilinear map). One can thus apply eigenexpansions 
etc. to, say, optical cavities j^. In particular, the di- 
agonal bilinear map, known to Zeldovich |2^| long ago 
at least for eigenf unctions, has been used for stationary 
perturbation theory (0 and references therein). 

This previous work points to a remaining subtlety. 
Namely, while completeness is a simple consequence of 
linear algebra and counting in the discrete case, the 
modes of the open cavity are complete only if V(x) or 
one of its derivatives (depending on the permissible de- 
gree of regularization in the field expansion) are discon- 
tinuous at x = a — demarcating the cavity. One should 
thus study how the eigenexpansion behaves as A — > 0, de- 
pending on whether or not this discontinuity condition is 
satisfied in the limit [one can speculate that some modes 
of (|3.8|) - (|3.9|) run off to —ioo if it is violated]. Ideally 
this should yield a completeness proof alternative to the 
existing one, based on the analyticity of the continuum 
Green's function for Imu < but limited to 1-d space. 



C. Wronskian 

Assuming nearest-neighbor couplings has turned these 
continuum models into second-order differential equa- 
tions, for which there is an extra ingredient. One can 
define, for any lo, two solutions f{x,u) and g(x, uj) sat- 
isfying the left and the right boundary condition respec- 
tively. For, say, a string fixed at both ends, /(w,0) = 0, 
g{uj, a) — 0. The x-independent Wronskian is J(u>) — 
g(x,uj)d x f(x,uj) — f(x,u))d x g(x,ui). At a zero u>j of J, 
f on g satisfy both boundary conditions, hence are eigen- 
functions. Thus, the Wronskian plays the role of the 
characteristic polynomial. Moreover, the (normalization- 
independent) Green's function can be represented as 
G(x<y, w) = f(x, u)g(y, lo)/J{uj) [G(y, x) = G(x, y)] . 

There are however differences from the general finite- 
dimensional case: 1-d second-order differential equations 
do not admit level crossing, so a multiple zero of J always 
implies criticality. Also, the one-sided function f (or g) 



VI. DISCUSSION 

In the conservative classical system's eigenvec- 

tors (e.g., electromagnetic plane waves) are a starting 
point for discussing nonlinear phenomena (e.g., second- 
harmonic generation), quantization (e.g., plane- wave cre- 
ation and annihilation, Feynman propagators), as well as 
quantum interactions (e.g., a vertex in a Feynman dia- 
gram linking several propagators). Similarly, our eigen- 
expansion for ohmically damped linear classical oscilla- 
tors is relevant for nonlinear and/or quantum systems 
as well — especially useful for systems with many degrees 
of freedom dominated by a few modes (typically at low 
temperatures or when some masses are large). Here we 
only provide a sketch; details will be given elsewhere. 



Nonlinear oscillators 



Consider for example a set of oscillators described by 

d?0(a) + [T(a,p)d t +K(a,f3)+\(a,0, 7 )ct>(j)}ct>(f3) = . 

(6.1) 

Using the expansion <p(t) = a ^ifi)fji one shows that 

{d t + Wj)a j (t) = -i\ jk i a k (t)a l (t) , (6.2) 
\ 3 u = A(a,/3, 7 )/ J (a)/ fc (/3)/ i (7) ■ (6.3) 

For small A, these can be solved perturbatively. In partic- 
ular, (|2.12|) allows a J (0) to be found by projecting 0(0). 



B. Quantum and thermal correlation functions 

To deal with quantum and/or thermal physics (for the 
moment ignoring nonlinearities), we need to modify two 
ingredients. First, the equation of motion becomes 

al^(a) + r(o, /8)dt0(/9) + K{a, /3)0(/3) = rj(o, t) . (6.4) 

The noise r\ coming from the bath satisfies the 
fluctuation-dissipation theorem plj 



2w 



1 - e-^/ T 



T(a,0), (6.5) 



where (• • • ) denotes the expectation value at tempera- 
ture T, and h = /cb = 1. For finite-T classical dynamics, 
the rhs reduces to 2TT(a, f3). In two-component form, 



(id t - U)4> = iS(t) , S = 



(6.6) 



can be studied near uij, of use for critical points 



ton / (or 
ints [2J. 



In the eigenvector basis [26j, 

(d t + = Si(t) ee (£, S(t)) = (f j \S(t)) . (6.7) 

Correlators such as (a^ (t)a k ) are then related to 
(S : >(t)S k ), which can be obtained from l|6.5|) and l|6.6|l . 
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Classically, (|6.4fl - H6.7fl can already yield nontrivial 
Langevin dynamics. In quantum mechanics, <fi and hence 
its coefficients a? become operators. Commutators are 
evaluated by first substituting a? = (/.,■,</>) [implying 
a~ J = ^fi(a^y when combined with 12.14fl ]. and then im- 
posing the canonical [4>(a), <j){f3)\ = iS(a,f3). One finds 



(a k )t] = ( Wj +«iJ)/ J -(a)/ fc *(a) = A^ 



(6.8) 



For r —> 0, the orthogonality under (|2.9[1 is (uij+Uk) x 
fj(a)fk(a) — <P , with all quantities real. Hence, A jfc = 
5 3 +0(T) so that a 3 and (aPy generalize the usual anni- 
hilation and creation operators. Crucially, however, their 
commutators l|6.8l) are not diagonal — rooted in the very 
nature of the bilinear map, and implying that the Feyn- 
man propagator is likewise non-diagonal. Diagrammatic 
expansions will thus show interesting features. 
In fact, we can calculate the correlation function 



(6.9) 



in a way that bypasses some of these technical complica- 
tions. For lo 6 R, from general considerations [25| 

F(a, (3, u) = 2(1 - e^^f 1 Im G QP (a, 0, w) . (6.10) 
Using l|2.24|l for Q then gives 

2z ^ufjWfM) 



1 



(6.11) 



with poles at both the uij and the Matsubara frequencies 
uj m = 2rmriT. Substituting the integrated 16.7fl into the 
eigenexpansion of JfUHJ also yields ffTTty . A gain, we omit 
details, and only stress that the eigenbasis provides a 
direct correspondence with the frequency singularities. 

The Feynman propagator can be evaluated in terms 
of F and Q . The perturbative treatment of interactions 
of Section IVI Al (readily promoted to the operator level) 
then opens the way to a diagrammatic analysis. 



C. Conclusion 

We have established an eigenvector expansion for a sys- 
tem of ohmically damped oscillators such that almost all 
the usual results for conservative systems carry over. In 
a sense, introducing the bilinear map is equivalent to in- 
verting the (potentially huge) matrix (fj\fk)- Initial- 
value problems are readily solved, and RSPT is valid 
for the complex eigenvalues. The possible extension to 
nonlinear and quantum phenomena has been sketched. 
The present formalism fails when any mode is critically 
damped; this is handled in the sequel p. 

These developments should be useful whenever there 
is (a) dissipation and (b) many degrees of freedom. Ex- 
amples would include mechanical vibrations, acoustics, 
electromagnetism in absorbing (as well as gain) media, 



large Josephson-junction arrays and macromolecule vi- 
brations in a viscous medium. 

An extension to non-ohmic damping (which, by the 
Kramers-Kronig relations, typically would entail disper- 
sion in the force constants as well) would enable modeling 
the many dispersive systems encountered in applications. 
Moreover, a closer study of the quantum dynamics re- 
quires a high-frequency cutoff on the damping. 
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APPENDIX A: FAST MODES 

1. Equations of motion 

If one particle is much heavier than all others, the lat- 
ter can be treated as "fast modes", the momentum of 
which is eliminated. In preparation for such a study and 
for another application in Appendix IA 31 here we consider 
one fast mode a = N+l of mass e <C 1. Specifically, 

d 2 t (j){a) + [T(a,(3)d t +K(a,P)]ct)(l3) = -A(a)<p(N+l) 

(Al) 

for 1 < a < N (also the range of the summation con- 
vention), while the light particle is coupled to the others 
only by force constants A(a), not by cross-damping |27j : 



(ed? + jd t + k)4>(N+1) = -A(a)(p(a) 



(A2) 



This model (easily extended to more light masses) is con- 
tained in (|1.1[) . if one rescales <f>(N+l) = y/e<f>(N+l). 

In the adiabatic limit e — > 0, the last particle's momen- 
tum follows the other variables: 

d t (f){N+l) = -af>{N+l) ~ B(a)4>{a) , (A3) 

with c = k/j and 

B(a) = -f- 1 A(a) . (A4) 

The dynamical system defined by I|A1[) . (|A3() . and (|A4() 

is dissipative even if T = 0: l|A3|) breaks invariance under 
t i—* —t. The system is stable if the original one is, i.e., 
if 



(A5) 
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Assuming K > 0, Sylvester's theorem 01 guarantees 
stability if det/C > 0; in particular, one needs C7 > 0. 
Since 7 is taken positive in (|A2|> one finds that c > 0, 
which is plausible given the form of (|A3|) . 

On the Hamiltonian level, (|A3I) means that (f>(N+l) is 
not dynamical but obeys a constraint instead, hence the 
phase space has 2N+1 dimensions. Thus, define 

= (0(1), . . . , <j>(N); 0(1), . . . , $(N); 0(A +1)) T . (A6) 

Then the dynamics takes the standard form l|2.2(l . with 

/ I 0\ 
W . = * -if -T -A . (A7) 

V-£? t - c y 

2. Bilinear map 

Since the above is merely a limiting case of our general 
model, the bilinear map should be 

(0, ip) = i [(t>(a)ip(a) + 4>(a)ip{a) 

+ 4(N+1)]j(N+1) + £(N+l)-ip(N+l) 

+ <p(a)T{a, (3)4>{(3) + $(N+l)e- 1 j<ip(N+l)] 

= i[4>{a)ijj{a) + 4>(a)ip(a) + cj>(a)T(a, P)tp(P) 
+ 0(/v+l) 7 <0(iV+l)] + 0(e) (A8) 

in terms of the vectors (|A6|) . By l|2.15fl . this leads to 

/r 7 0\ { K A\ 

g.. = i [I , = -7 . (A9) 

\0 7/ \,4 T 07/ 



Interestingly, when starting out by postulating (|A3I) in 
addition to QAl[l . the proportionality (|A4() follows from 
demanding that l|A9|l be symmetric. Likewise, requir- 
ing the existence of a non-increasing energy function 
bounded from below then leads to l|A5(l . 

The adiabatic limit can similarly be taken in essentially 
all results of the previous sections, yielding a realization 
of our framework in an odd-dimensional phase space. 



3. Example 

Because of <|A4fl , we should demonstrate nontrivial ex- 
amples. Consider a string of particles < a < 7V+1, 
spaced by A, with a = (N+l)A and 0(0) = 0. Let the 
force constants be given by (|3.4[) with fcA 2 = 1, while for 
simplicity F = 0. The nearest-neighbor couplings l|3.4() 
also act between particles and 1, and between particles 
N and 7V+1. The latter implies A(a) = -A~ 2 5(a, N). 

The right boundary condition is the discrete version 
of IJ53J, namely d t <j>(N+\) = -A" 1 [<p{N+l) - 0(A)]. 
Compared to (|A3|) . we have B(a) = —A~ 1 S(a, N) and 
c = A -1 . Crucially, l|A4|l is satisfied, with 7 = A -1 . 

At least formally, for N — > 00 (with a fixed) one re- 
trieves the model of Section IV Bl Indeed, the continuum 
limit of l|A8(l agrees (up to an overall factor of A) with 
(|5.()|) . It would be instructive to compare the convergence 
of the present discrete model with the more physically 
motivated one from the main text. 
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